Ion-induced nucleation in polar one-component fluids 
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We present a Ginzburg-Landau theory of ion-induced nucleation in a gas phase of polar one- 
component fluids, where a liquid droplet grows with an ion at its center. By calculating the density 
profile around an ion, we show that the solvation free energy is larger in gas than in liquid at the 
same temperature on the coexistence curve. This difference much reduces the nucleation barrier in 
a metastable gas. 



I. INTRODUCTION 



Ion-induced nucleation in water vapor plays a decisive 
role in many atmospheric phenomena. It is well-known 
that the nucleation rate is much enhanced in polar flu- 
ids in the presence of ions^'^. However, precise experi- 
ments on this effect have been difficult despite the long 
history of this problem since Wilson's cloud chamber 



experiment 
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Theoretically, we need to understand 



statistical behavior of polar molecules around an ion^ and 
heterogeneous nucleation triggered by such ion-dipole in- 
teraction in metastable states. This problem is thus of 
fundamental importance in physics and chemistry, but it 
has rarely been studied in the literature. 

Originally, Thomson^ calculated the chemical poten- 
tial of vapor in the vicinity of a charged particle using 
continuum electrostatics and predicted formation of a 
small liquid droplet with the ion at its center. How- 
ever, its radius is only SA in room-temperature water 
vapor, so Thomson's picture is an oversimplified one. Re- 
cently, ion-induced nucleation has been studied using a 
density functional theory^'^, a Monte Carlo method'*, and 
molecular dynamics simulations^ . These numerical stud- 
ies have shown clustering of polar molecules around an 
ion and alignment of the dipoles along the electric field. 
After the statistical average, the first effect gives rise to 
a density inhomogeneity n{r) or electrostriction around 
an ion. It is worth noting that Born calculated the solva- 
tion free energy due to an ion in continuum electrostat- 
ics, accounting for polarization alignment (in the linear 
response) but neglecting electrostriction''''^^'^^. 

The aim of this work is to present a Ginzburg-Landau 
theory of ion-induced nucleation in polar one-component 
fluids in a gas phase. It will be based on our recent 
theory of solvation effects in near-critical polar binary 
mixtures^^'^'^. Our theory will take account of both elec- 
trostriction and polarization alignment. As will be dis- 
cussed, the Born formula becomes a poor approximation 
in strong electrostriction, particularly when an ion is in a 
gaseous polar fluids. Merits of our approach are its sim- 
plicity and its ability to describe mesoscopic effects such 
as nucleation. In Section 2 a theoretical background of 
solvation will be presented. In Section 3 the free energy 
increase to create a critical droplet around an ion will be 
calculated in a gaseous polar fluid. 



II. THEORETICAL BACKGROUND 

A. Electrostatics 

We place a charged particle (ion) in a polar one- 
component fluid. We assume that the static dielectric 
constant e = e{n) is an increasing function of the fluid 
density n at each temperature T . In the continuum ap- 
proximation the electric field E = — V$ is induced by 
the electric charge and the electric potential $ satisfies 



V • eV* = -47r/9(r), 



(2.1) 



where p{r) is the charge density. Here the polarization 
P is assumed to be aligned along E as P = with 
X = (47r)^^(e — 1) being the polarizability^^. We take 
the origin of the reference frame at the center of the ion. 
It is convenient to assume that the charge density is ho- 
mogeneous within a sphere r < i?i as 



p = Ze/iAnRf/i) {r < i?;) 



(2.2) 



and vanishes outside the sphere r > Ri. The total charge 
is given by Ze. In our theory Ri is a phenomenological 
parameter representing the ion radiuS"'^'*. Then n = n{r) 
and <& — <I>(r) depend only on r = \r\. In this case 
Eq.(2.1) is solved to give 



An 



drirjpin) = ^e{r), (2.3) 



where ^'(r) = d'^{r)/dr, 0{r) 1 for r > Ri, and 

e{r)^{r/Rif (r<i?i). (2.4) 
The electrostatic free energy is simply of the form'^-'^. 
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K = / dr—eE^ = -Z^e^ 



dr—0{ry. (2.5) 



Its functional derivative with respect to n yields the 
chemical potential contribution arising from the electro- 
static interaction, 



on 



ZW_9{r 

87re2 r' 



(2.6) 



where s' = {de/dn)T- This quantity is negative and 
grows strongly as r approaches Ri, leading to accumu- 
lation of the fiuid particles around the ion. Here we 
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are neglecting the nonlinear electric field effect near the 
ion^'^^'^^, for which see a discussion in the summary. The 
effect arising from the complex molecular structure of 
solvent molecules is also beyond the scope of this work, 
which leads to the dependence of the nucleation rate on 
the sign of the charge of ions^ . 

A number of theories have been presented to describe 
the overall density and temperature dependence of the 
dielectric constant e = e{n,T). Recently, Fernandez et 
al}"^ analyzed a wide range of data of e for water and 
steam using the formula proposed by Harris and Alder^^ 
for polarizable polar fluids, 



£-1 

7+2 



= A-nn 



OL-n 



+ 



kBT (2£ +!)(£ + 2) 



(2.7) 



where a„i represents the molecular polarizability and /x is 
the dipole moment of a molecule(~ 2D for water). This 
form reduces to the Clausius-Mossotti formula in the ab- 
sence of a permanent dipole and tends to the Kirkwood 
formula (e — 1)(2£ + !)/£ = Anng^^ /k^T for nonpolariz- 
able molecules with permanent dipole momc nt^''. Here g 
is the so-called Kirkwood correlation factor arising from 
the correlation between the molecular orientations due 
to nondipolar interactions (2—3 for liquid water), so g 
weakly depends on n as revealed in experiments^". In 
Fig.l we plot £ as a function of n at T/Tc = 0.675 
(T = 437 K) and 0.928 (T = 600 K) below on the basis 
of the formula (2.7) using data of am and gjj,^ for water 
in Ref.l7. In our numerical analysis the quantity £'/£^ 
in Eq.(2.6) will be calculated using the formula (2.7). 



B. Ginzburg-Landau free energy 

For simplicity we use the Helmholtz free energy density 
/ = fi^jT) in the van der Waals theory^^. 



/ = fceTnln 



1 — VqU 



— k-QTn — Ac^^vqu 



(2.5 



where and Cvw represent the molecular hardcore vol- 
ume and the strength of the attractive part of the pair 
interaction, respectively, and Ath = h{27r/m,kBTy^'^ is 
the thermal de Broglie length. It follows the well-known 
expression for the van der Waals pressure, 

p = nn-f=- 4evwVon^, (2.9) 

1 — VqU 

where /U = df /dn is the chemical potential per particle 
of the form. 



IJL= fcBTln 



1 — VQtl 



+ 



kBT Von 



^von. (2.10) 



1 — Von 

The isothermal compressibility Kt behaves as 

nkBKT = {l-vonf/[T-T,{n)], (2.11) 



where Ts{n) is the spinodal temperature given by 

T,(n) = 8evwWon(l - vonf/kB (2.12) 

Maximization of Ts{n) with respect to n gives the critical 
temperature and density. 



Te = 32evw/27fcB, ric = l/Swo- 



(2.13) 



In Fig. 2 we show the phase diagram in the T-cj) plane 

for the van dor Waals fluid, where <f> = vqu is the nor- 
malized density or the volume fraction of the hard-core 
region. In this work we will set vq = 30.96 x 10~^^ cm^ 
and Cvw/fce = 545.98 K, which follow from the critical 
values, Uc = 1.076 x lO^^cm-^ and = 647.096 K, for 
water^^'^^. We introduce the van der Waals radius. 



1/3 



(2.14) 



which is 3.1 A for water. 

We assume that n tends to a homogeneous value n^o 
far from the ion r ^ oo and that the temperature T is 
homogeneous throughout the system. These assiimptions 
are allowable in calculating the equilibrium solvation pro- 
file and the critical droplet in nucleation^^. The chemical 
potential and the pressure far from the ion are written 
and Poo , respectively. We consider the grand po- 
tential f2 (equal to —pV for homogeneous states) at the 
temperature T and the chemical potential fi = /^oo- Its 
increase due to an ion including the gradient and electro- 
static contributions is written as 



Ail 



dr 



(215) 



where C is a constant and 

g{n) = /(n) - iJaon + Poo 

= /("-) - /oo - ^^oa{n - rioo), 



(2.16) 



with /oo being the Helmholtz free energy density far from 
the ion. The AO will be called the solvation free energy. 
Here g and dg/dn tend to zero at n ^ so that we have 
g = {n — noo)^/2n1^KT{noo) far from the ion. It is worth 
noting that if the electrostriction is neglected or n = rioo 
in the whole space, AO consists of the electrostatic free 
energy Fg only, leading to the Born contribution^i. 



AORorn = dr 



Ze 
Sir V £00?"^ 



2^2 



2£oo-Rb 



(2.17) 



where £00 = £("oo)- Here Rb = Ri if the space integral is 
outside the ion (r > Ri)^^ with the upper bound pushed 
to infinity. However, this contribution much overesti- 
mates Fe when n(r) approaches a liquid density close to 
the ion and tends to a smaller (gas or supercritical) den- 
sity far from it. In such cases F^ is of order Z'^e^ /2e(RB 
with ££ being the dielectric constant in liquid^'^, which is 
much smaller than AOBom in Eq.(2.17) for £oo- 
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In Eq.(2.15) Fe arises from the ion-dipolc interaction 
among the ion and the solvent molecules. However, there 
is also a repulsive interaction among them. To account 
for it we may include a potential Vis (r) to obtain another 
form of the grand potential increase given by 

Afi' = Aa + j drvis{r)n{r), (2.18) 

as in the density functional theory in Ref.6. The poten- 
tial Wis(r) should grow strongly close to the ion center. In 
this paper we use the truncated Lenard- Jones potential 
Vis defined by 

Vis{r) = evw[(tTis/r)i' - {ai^/rf] (r < a;,), (2.19) 

and Vis{r) = for r > ctis with = a/2 + Ri. We will 
show numerical results starting with either Eq.(2.15) or 
Eq.(2.18). 

In this paper the coefBcient C of the gradient term in 
Eq.(2.15) will be set equal to 

C = lOo-^fceT. (2.20) 

As will be shown in Appendix A, this choice of C is 
consistent with a calculation of the surface tension 7 
for water in Ref.23. The thermal correlation length 
^ = {Cu'^KtY/'^ is then about a for T/T^ = 0.675 and 
about 2.2(7 for T/Tc = 0.928 in liquid on the coexistence 
curve. 



In numerically solving Eq.(2.21) the ion radius will be 
assumed to be given by 

Ri = 0.3a. (2.24) 

which is taken to be considerably smaller than a and 
then (Tis = 0.8(7 in Eq.(2.19). In the litciraturc'"''!'' the 
ion radius has been estimated for various ions in water. 
For example, in Table V in Ref.5, the bare ion radius is 
approximately 0.68A for Li+, 0.95A for Na+, 1.69A for 
Cs^+, and 0.5 A for Al^+, which are considerably smaller 
than the corresponding hydration shell radius. The bare 
radius in Ref.5 corresponds to Ri in our theory. 

In Fig. 3 we show the normalized density profile (f>{r) = 
Von{r) around an ion in the gas phase case (f>oo — v^noo = 
0.049 for A = 5 and 50 by minimizing Ml' in Eq.(2.18) 
including Vis- Note that n(r) tends to zero rapidly at 
small r. Similar density profiles were calculated in the 
previous studies^'^'®. On the other hand, if we minimize 
Afi in Eq.(2.15) without Vi^, n{r) tends to a positive con- 
stant n(0) at small r < as for the concentration in 
our previous work^'^. We show that the solvation den- 
sity profiles change continuously close to the coexistence 
curve with a change in the temperature difference, 

6T = T- Tex, (2.25) 



C. Solvation in gas 

When an ion solvates the surrounding fluid in stable or 
metastable one-phase states, we may determine the den- 
sity profile n = n{r) by requiring the extremum condition 
SAQ/6n = or SAQ'/6n = expressed as 

fi-CV^n- iJ,oo = €yy,Aa^, (2.21) 

where jl is equal to /x in Eq.(2.10) for 6Afl/dn = or to 
/i + Vis for SAfl' /Sn = 0. The right hand side is /ig in 
Eq.(2.6). We introduce the dimensionless strength of the 
electric field, 

Z2e2 4 Z^e^ 
87r(7evw 277r(7fcBTc' ^ ' ^ 

using in Eq.(2.8) or T,. in Eq.(2.13). Here we es- 
timate A — 3.9Z^ for water. As shown in our previous 
work^'^, the gradient term becomes negligible in Eq.(2.21) 
far from the ion. The long distance behavior of the devi- 
ation 5n = n — rioo far from the critical point is 

8n ^ D^e^^Aa/r'^ (r > a) (2.23) 

with being the value of r^Kq-e' je^ far from the ion. 
The solvation profile is more complicated near the critical 
point due to the growing correlation length^^. 



where Tex is the coexistence temperature at given density. 
This means that there is no discontinuous change in the 
solvation profile and the free energy near the coexistence 

curve^®. 

Let us consider the condition of strong electrostriction 

in the gas phase, where n{r) assiimcs a liquid density 
close to an ion. Note that the chemical potential on the 
left hand side of Eq.(2.21) is of order k^Tj (1 — v^n) and 
e' is of order v^e in liquid states with fceT ~ e^,^ and 
n Vq^ . Thus strong solvation condition is given by 

A(<7/iii)V£^ > 1, (2.26) 

where £^ is the dielectric constant in liquid. The degree 

of solvation strongly increases with decreasing the rel- 
ative ion size Ri/o, as is well-known in the literature^. 
The shell radius RsheW may be defined such that n(r) is 
on the order of the liquid density within the shell region 
r < i?sheu- Then we find A{a/Rsheii)'^/£e ~ 1 to obtain 

Rsheii ~ (7{A/sey^\ (2.27) 
which can exceed a for A/ee >1. 
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D. Solvation on the coexistence curve 

Due to the large density variation around an ion in gas, 
the effect of electrostriction is more marked in gas than 
in hquid. As a result, the solvation free energy is larger 
in gas than in liquid on the coexistence curve. We show 
/3Ail without I'is in Fig. 4 and /3Ail' with in Fig. 5 for 
the two phases along the coexistence curve. We notice 
that /3Af2' even exceeds 1000 on the right panel of Fig. 5. 
This is because e{n) decreases from a liquid value to a 
gas value in the region r < ais in the presence of Vis- As a 
result, the integration of the electrostatic energy density 
in the region Ri < r < (Tig becomes of order AilBorn in 
Eq.(2.17) (~ lO^fcBT for A = 50). 

In the next section we shall see that relevant in nu- 
cleation is the difference between the gas value and the 
liquid value of the solvation free energy, 

Asol = {An)n=n^ - (A0)„=„, without Vis 

= (Ar!')n=n, - (A0')„=„, with Vis. (2.28) 

Here the grand potential increases, Af2 and Af2', are cal- 
culated using the inhomogeneous solutions of Eq.(2.21), 

n = ng{r) for gas or n = nf {r) for liquid on the coexis- 
tence curve. For the case with Vi^ we define the grand- 
potential density, 

uj'{r)=g + Vi, + ^\Vn\'' + ^E\ (2.29) 

Then we have 

/•OO 

A,oi = 477 / drr^iu'Jr) - u'^{r)], (2.30) 

where n = ng(r) in lo'^ and n = ne{r) in w^. The two 

terms, i-^gir) and (-o'g{r), in the integrand mostly cancel 
in the region r < (Tig. In Fig. 6 we show PAgoi on the 
coexistence curve for A = 5 (left) and 50 (right). As ex- 
plained below Eq.(2.17), the Born approximation much 
overestimates Agoi in strong solvation. 



III. ION-INDUCED NUCLEATION IN GAS 

A. Critical droplet and KO approximation 

In nucleation theory^^'^^'^* crucial is the free energy in- 
crease to create a critical droplet, which will be written 
as Wc and will be called the nucleation barrier. Its def- 
inition for ion-induced nucleation will follow in Eq.(3.3) 
below. The nucleation rate J is the birth rate of droplets 
of the new phase with radius larger than the critical ra- 
dius i?c emerging in a metastable state per unit volume. 
It is approximately given by 

J= Joexp(-/3We), (3.1) 



independently of the details of the dynamics. Here Jo 
is a microscopic coefficient (of order n/rm with Tm be- 
ing a microscopic time away from the critical point) and 
(3 = l/ksT. A unique aspect of nucleation in gaseous 
polar fluids is that the solvated region around an ion 
can serve as a seed of a nucleating liquid droplet in a 
metastable state. This leads to a considerable decrease 
in Wc and a dramatic increase in J in the presence of a 
small amount of ions. 

Mathematically, the extremum condition (2.21) holds 
also for the critical droplet profile, 

n = ncri(r), (3.2) 

which assumes liquid values for r <Rc and tends to a 
metastable gas density Uoo at large r^Rc. This solution 
of Eq. (2.21) is unstable; that is, if we superimpose a small 
density perturbation inside a critical droplet, AQ (or 
An') decreases^^. In Fig. 7 with Vis being included, the 
normalized density 4){r) = vancri{r) for a critical droplet 
is displayed around an ion with A = 50 for T/T^ = 0.928 
(left) and 0.675 (right). The dotted and broken lines rep- 
resent the solvation profile in a metastable gas and that 
in a stable hquid, respectively. Here dT/T^ is —0.011 
(left) and —0.046 (right), so Rc is larger for the former 
case. 

The nucleation barrier Wc is then written as 

Wc = (Af2)„=„„j - (An)n=n^ without Wis 

= (A0')«=«„. - (Ar!')„=«, with Vis, (3.3) 

where ncri{r) is the critical-droplet solution and ri.g(r) is 
the inhomogeneous solution for the metastable gas. In 
terms of w' in Eq.(2.29) we have 

/•oo 

Wc^Att drr^luUr) " ^^(O], (3-4) 
Jo 

for the case with Vis- Here n = ncii(r) in cj^j.j and 
n = ng{r) in Wg. The numerical calculation of Wc in 
this manner will be referred to as KO approximation. 
In Fig.7 the calculated f3Wc is 50.58 (left) and 101.24 
(right). In Fig. 8 we show the dimensionless quantity 
4:TT(3ar^[uj'^^.{r) - ^(r)] for the case of T/Tc = 0.928 
(left) in Fig.7, which is the integrand in Eq.(3.4) mul- 
tiplied by f3a. It consists of a negative part within the 
droplet r <Rc and a positive surface part at r ~ i?c, and 
its integral with respect to r/a is equal to pWc ~ 50. 

B. KO' approximation 

In the KO approximation for Wc we need to calculate 
'^cri('') for each given set of the parameters ST, A, and Ri. 
We here propose a simpler (KO') approximation valid for 
shallow quenching. That is, we set 

Wc = -Asoi + W°. (3.5) 
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where Agoi is the constant independent of R defined by 
Eq.(2.28) and VK" is the nucleation free energy in the 
absence of ions or for homogeneous nucleation. The first 
term arises from the difference between the solvation free 
energy in gas and liquid, which may be equated to the 
values on the coexistence curve for shallow quenching as 
indicated by Fig. 3. In curves in Figs. 9 and 10 to follow, 
W° will be numerically obtained from Eq.(2.21) for pure 
fluids or for Ze = and vis = 0. 

For large droplets {R 3> £,) we may propose the follow- 
ing droplet free energy, 

Wko{R) = -Asoi + 47r7i?2 _ —dijR\ (3.6) 

The last two terms have the meaning of the minimum 
work needed to form a droplet with radius R in homo- 
geneous nucleation^^. Namely, the second term is the 
surface part with 7 being the surface tension. The third 
term with S/j, = —g(ne) is the bulk part, which is negative 
or > in metastable states. Notice that the dotted 
curve in Fig. 8 represents — 47r(/3fT5/i)r^. See Appendix B 
for more discussions on 5^,. This Wko(-R) is maximized 
at the Kelvin radius, 

i?K = 27/(5/x, (3.7) 

which is the critical radius in homogeneous nucleation. 
The maximum value is the nucleation barrier. 

We ^ -A,oi + ^7r7i?2 . (3.8) 

In this simple case Eq.(3.1) may be rewritten as 

5^l - (16^/37V3)i/V[ln(./o/ J) + dKoiY'^ ■ (3.9) 

Appreciable droplet formation is observed for sufficiently 
large nucleation rate J dependent on the experimental 
method (~ 1/cm^sec, say) and the above formula deter- 
mines the so-called cloud point of nucleation^ ^. 



where the upper bound of the integration is R since e(n) 
changes from eg to at r = i? with increasing r. The 
first term in Eq.(3.10) arises from the lower bound and 
turns out to be the Born approximation of — Agoi in 
Eq.(2.28). 

Thomson^ found that the sum of the second and third 
terms on the right hand side of Eq.(3.10) is minimized at 
the Rayleigh radius"^^, 

Rt^ = {a/Sn^y/^. (3.13) 

He concluded that an ion in gas is surrounded by a spher- 
ical liquid region with radius However, as estimated 
by Thomson himself, i?R is only about SA for an ion with 
Z = 1 in room-temperature water vapor, while Eq.(3.10) 
should be valid only for large R. For not very large Z we 
can see that i?R cannot exceed even the interface thick- 
ness at any temperatures^^. In addition, use of the sur- 
face tension makes his argument applicable only in the 
region ST <0. As shown in Fig. 3, the solvation density 
profiles of small ions should be continuous through the 
coexistence curve. Obviously, the second term (oc Rr^) 
in Eq.(3.10) is negligible for R i?R and the Thom- 
son theory is not well justified (if applied to microscopic 
ions) . 

We explain how the nucleation free energy Wc is cal- 
culated in the Thomson theory^. The extremum con- 
dition dWT{R)/dR = becomes a;^(l - x) = a* with 
X = R/Rk, where 

a* = a/8n^Ri = {Rn/RKf. (3.14) 

We can see that Wt{R) has a maximum Wmax and a 
minimum Wmin for 

a* < 27/256 or Rr/Rk < 3/4''/^ (3.15) 



C. Thomson approximation 

Since Thomson's work^, the droplet free energy for 
large R has been assumed to be of the form^'"^"'"^^, 

Wt{R) = -^.+^ + 47r7i?' - Y^t^R", (3.10) 

where a is written in terms of the dielectric constant in 
gas Sg and that in liquid as 

a=^Z'e\l/e^-l/e,). (3.11) 

The first two terms in Eq.(3.10) represent the difference 
of the electrostatic free energy F^. between the two phases, 

(AFe)Born = - / dr^, (3.12) 
JRi 



Since i?R is microscopic for not large Z, wc have a* ^ 1 
for shallow quenching. Obviously, the radius giving the 
maximum is nothing but the critical radius Re- For small 
a*, Rc behaves as Rc = i?K(l — a* + ■ ■ while the ra- 
dius giving the minimum is nearly equal to _Rr. If Thom- 
son's idea is extended, a preexisting liquid droplet with 
radius i?R, grows into a larger liquid droplet in nucleation. 
Within his theory, the nucleation barrier is given by^ 

We = W„,ax - W^in 

= ^TT'yR^{l + 8a*)-12TT'yRl^ + ---. (3.16) 

In the first line the first term in Eq.(3.10) is canceled 
to vanish. The second line is an expansion in powers of 
{a*y^'^ for small a* . This approximation was adopted in 
Ref.3, but it is not well justified for small i?R and hence 
for small ions. 
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D. Numerical results for Wc 

Without Vis in Fig. 9 and with Vis in Fig. 10, the nor- 
niaUzcd nucleation barrier f3Wc is plotted as a function 
of A at T/Tc = 0.928 (left) and 0.675 (right). The pro- 
files of the critical droplets are shown in Fig, 7 at ^ = 50 
for these temperatures. The points (•) (KO) are ob- 
tained from the profile of the critical droplet satisfying 
Eq.(2.21). The solid line represents /3Wc in Eq.(3.5) 
(KO'), while the broken line represents pWc in Eq.(3.16) 
(Thomson). It is remarkable that there is no essential 
difference between Fig. 9 without Wis and Fig. 10 with Vis- 
(i) For T/Tc = 0.928, good agreement is obtained be- 
tween the KO result and the KO' result because of rel- 
atively large Rc, while the Thomson result is larger by 
about 10. For example, for T/T^ = 0.928 and A = 50 
in Fig.lO, (3Wc is given by 50.58 (KO), 50.71 (KO'), and 
61.59 (Thomson), (ii) For T/T^ = 0.675, the KO val- 
ues of agree with the KO' values for A <10, but 
become considerably larger for A >20. In fact, /3Wc is 
given by 101.24 (KO), 77.80 (KO^, and 80.96 (Thom- 
son) for T/Tc = 0.675. In this lower temperature case 
Rc is of order iJsheih as shown in Fig. 4, and hence the 
expression of Wt{R) in Eq.(3.10) should not be a good 
approximation. 

For these two temperatures /3Wc decreases with in- 
creasing A. At small A{<5) the curves are rather steep, 
mainly due to the small size of Ri/a = 0.3 in Eq.(2.24). 
In the density functional theory Kusaka et al.^ calcu- 
lated Wc for i?i = cr as a function of another parameter 
X (= Zepo/a^kBTc with po being the dipole moment of 
a polar molecule) and found that it is rather weakly de- 
pendent on X for small x- 

IV. SUMMARY AND CONCLUDING REMARKS 

In summary, wc have presented a simple theory of ion- 
induced nucleation in polar one-component fluids, where 
the effect of electrostriction is crucial. We summarize our 
main results. 

(i) We use a simple continuum model, where the elec- 
trostatic free energy is given by Eq.(2.5) and the grand 
potential increase due to an ion is given by Eq.(2.15) or 
Eq.(2.18). Particularly, we use the Harris- Alder formula 
(2.7) for the dielectric constant as a function of the den- 
sity and the van der Waals form for the Helmholtz free 
energy density in Eq.(2.8). The density profiles around 
an ion are numerically obtained as illustrated in Figs. 3 
and 7. 

(ii) We have proposed the droplet free energy Wko{R) 
in Eq.(3.6) for a large liquid droplet in a metastablc gas 
with an ion at its center. It is remarkable that the dif- 
ference of solvation in gas and in liquid gives rise to the 
negative background contribution — Agoi, which much fa- 
vors nucleation around an ion in a metastable gas. 

(iii) The nucleation barrier Wc has been calculated in 



three manners. In our Ginzburg-Landau scheme we have 
directly sought an unstable solution of Eq.(2.21) to ob- 
tain the points (•) (KO) in Figs. 9 and 10. It may also 
be calculated more approximately using Eq.(3.5) (KO') 
or the first line of Eq.(3.16) (Thomson). They yield the 
solid lines and the dotted lines in Figs. 9 and 10, respec- 
tively. We have found that the KO' result is close to the 
directly calculated result (KO) for relatively large critical 
radius Rc if the backgroimd is chosen to be — Agoi- 
(iv) We have presented the density profiles in Figs. 4 and 
7, which go to zero at the ion center owing to the pair 
potential Wig. However, we have obtained essentially the 
same results for the free energy barrier Wc in Fig. 9 with- 
out Vis and in Fig.lO with Vis- 
Finally we make some remarks. 

(i) We have used the Ginzburg-Landau theory with the 
simple gradient free energy. However, far from the crit- 
ical point as in the case of T/Tc = 0.675, the density 
functional theory^'^^'^^ with a nonlocal pair interaction 
would give more improved results. 

(ii) Because of its simplicity we have used the macro- 
scopic linear dielectric formula (2.7) even very close to 
the ion. This has been the approach in the previous 
continuum theories^'^", but it is not well justified. In 
particular, we should examine the effect of nonlinear 
dielectric saturation in solvation, which gives rise to a 
decrease in the effective dielectric constant close to the 

(iii) The difference Agoi of the solvation free energies in 
Eq.(2.28) is a relevant parameter generally in ionic solu- 
tions in two-phase states, although the present work has 
treated the case of a single ion. For example, it gives rise 
to ion-density differences inside and outside a wetting 
layer. 

(iv) Obviously Wc tends to zero on approaching the crit- 
ical point. However, this is the result for a single ion. 

We point out that even a small concentration of ions can 
drastically alter the phase behavior of near-critical polar 
fluids". 

(v) We can also construct a theory of ion-induced nucle- 
ation in binary fluid mixtures^^. 
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Appendix A 

Here we consider the surface tension 7 in the absence 
of ions in our model in the mean fleld theory^^. To have a 
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planar interface we assume that the fluid is on the coex- 
istence curve or the pressure is given hy p = Pcx{T) with 
the temperature below Tg. The grand-potential density 

g{n) = f{n) - ^ic^n + Pcx (A.l) 

is minimized at n = rig and rijt on the coexistence curve 
with the minimum value being zero. For our van der 
Waals free energy / in Eq.(2.8) and the gradient free 
energy in Eq.(2.15), a planar interface profile n = n{x) 
changing along the x axis satisfies dg/dn — Cd'^n/dx'^, 
yielding g = C{dn/dxY /2. Then 7 is expressed as 

/oo 
dx[g + C{dn/dxf/2] 
-00 

= (2C)i/2 / dn[g{n)]^/^. (A.2) 

In the vicinity of the critical point it follows the power 
law behavior, 

7 = aevw(l-T/T,)3/Va2, (A.3) 

where Cg is a constant. Widely in the range T/Tc >0.6, 
however, we have found that 7 numerically calculated 
from the second line of Eq.(A.2) with Eq.(2.13) can be 
fitted to Eq.(A.3) within a few percents if we set 

Cs = 3.0. (A.4) 

With this Cs Eq.(A.2) yields 7 ^ 235(1 - T/Tcf/"^ 
dyn/cm where we use the values of eyw and a given below 
Eq.(2.13) for water. 

Experimentally^'^, the surface tension of water behaves 
as 215(1 - T/Tcf" dyn/cm with v ^ 0.63 in the range 
1 — T/Tc <0.1. However, the calculated 7 fairly agrees 
with the experimental surface tension of water far from 
the critical point. For example, at T/Tc = 0.675, our 
7 is 42.5 dyn/cm, while the experimental value is 44.6 
dyn/cm. Since Cg oc C^/^, our choice of C in Eq.(2.21) 
is justified. 

Appendix B 

We examine the behavior of 5^1 = —g{ni) in Eq.(3.6) 
first fixing T below Tc- From Eq.(2.16) it becomes 

5ii = pt-Poo+ neiHoc - 
= -{fine) - Mcxn^ +Pcx) 

+ne{iJ.oo - Mcx) - (Poo - Pcx). (B.l) 

The quantities with the subscript £, 00, and cx arc those 
within the liquid droplet, in the metastable gas, and 
on the coexistence curve, respectively. See Ref.28 for 
a derivation of the first line of Eq.(B.l) from statistical- 
mechanical principles. In the second line, the first term 
is negligible since ne is close to the liquid density on the 



coexistence curve. In the limit Pcx —Poo ^ 0, the Gibbs- 
Duhem relation gives 

^ (An/ng)(poo - Pcx), (B.2) 

where An = — Ug. This relation is useful near the 
critical point. For dilute gas we have p = ksTn and 
Moo - Mcx = dpn{p)-'^ ^ kBTln{poc/Pcx) at fixed T. 
If rig <^ rii, the second term is much smaller than the 
first on the right hand side of Eq.(3.10), leading to the 
well-known expression^'^, 

fcBTln(poo/pcx). (B.3) 

In real experiments the temperature T^o in gas (written 
as T in the above discussion) may also be changed. Let 
us consider a reference state at temperature Tq and pres- 
sure po = Pcx{To) on the coexistence curve. In such cases 
we should replace Poc — Pcx(2o) in the above formulae 
by Poo - Pcx{To) - {dp/dT)cx{Too - To) =Poo- Pcx(T'oo) 
for shallow quenching, where {dp/dT)cx is the deriva- 
tive along the coexistence curve^^. In Fig. 2, for ex- 
ample, the temperature is changed across the coexis- 
tence curve at fixed average density no or at fixed sys- 
tem volume, where we may set Tq — Tex (no). This 
means Poo - Pcx(To) = {dp/dT)n{Too - To) in the 
gas, leading to 5^ = [Anidn/dT),x/nlKT]{To - Too) 
in terms of the compressibility Kt in the gas. Here 
use is made of the thermodynamic relation [dp/dT)n — 
{dp/dT)cx = —{dn/dT)cx/nKT on the coexistence curve 
(see Eq.(2.2.39) of Ref.[21]). 
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FIGURE CAPTIONS 



FIG. 1. Static dielectric constant e as a function of (p — vqu at T/Tc = 0.675 and 0.928 on the basis of the Harris- Alder 
formula (2.7) using data of water in Ref.l7. The dotted part corresponds to the two-phase region in the van der Waals 
model and is given by 0.0358 < (j)< 0.729 for T/T^ = 0.675 and by 0.168 < (j)< 0.519 for T/Tc = 0.928. 

FIG. 2. Phase diagram in the plane of T/Tc and cj) = Von with coexistence curve (solid line) and spinodal curve 
(dotted line). Use is made of the van der Waals model. At the points (•) numerical calculations are performed in this 
work. 

FIG. 3. Normalized density (p{r) = von{r) around an ion in the gas phase (t^oo = voHoo = 0.049) close to the 

coexistence curve (Tex = 0.721) for A = 5 and 50. The data arc for the three points with (f) = 0.049 in Fig. 2. The 
curves for T/Tc = 0.759 represent stable profiles, while the others metastable ones. However, for T/Tc = 0.591 and 
A = 50, there is no metastable solution with limr_>(x) 'n{r) = Uoo- 

FIG. 4. Normalized solvation free energy (3AQ without Vis around an ion for gas and liquid on the coexistence curve 
for A = 5 (left) and 50 (right). 

FIG. 5. Normalized solvation free energy 0Afl' with Vig around an ion for gas and liquid on the coexistence curve for 
A = 5 (left) and 50 (right). 

FIG. 6. Normalized solvation free energy difference (3Asoi defined by Eq.(3.5) along the coexistence curve with and 
without Vis for ^ = 5 (left) and 50 (right). 

FIG. 7. Normalized density (j){r) around an ion with ^ = 50 in a metastable gas with (poo = 0.18 (dotted line) and 
in the corresponding liquid (broken line) at T/Tc = 0.928 (left) and 0.675 (right). Also shown is the critical droplet 
obtained from Eq.(2.19). Here Vis is present and n(r) tends to zero at small r. 

FIG. 8. 4:'K(3ar^[uJc^i{r) — Wg(r)], whose integration with respect to r/a gives (iWc as in Eq.(3.4). The dotted line 
represents —4:TT{(3a6iJ.)r^ where appears in Eq.(3.6). The difference between these curves within the droplet is the 
solvation contribution. 

FIG. 9. Normalized free energy 0Wc without Vi^ at T/Tc = 0.928 (left) and 0.675 (right) as a function of A. The 
points (•) (KO) are obtained by solving Eq.(2.21). The solid line (KO') is calculated from Eq.(3.14), while the broken 

line (Thomson) is from the first line of Eq.(3.18). 

FIG. 10. Normalized free energy /?Wc with v.s T/T,, = 0.928 (left) and 0.675 (right) as a ftmction of A. The points 
(•) (KO) are obtained by solving Eq.(2.21). The soUd Unc (KO') is calculated from Eq.(3.14), while the broken hne 
(Thomson) is from the first line of Eq.(3.20). 
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